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Introduction 

Most practical combustion processes, as well as fires and explosions, exhibit some 
characteristics of turbulent diffusion flames. For hydrocarbon fuels, the presence of 
soot particles significantly increases the level of radiative heat transfer from flames. 
In some cases, flame radiation can reach up to 75% of the heat release by combustion 
(ref. 1). Laminar diffusion flame results (ref. 2) show that radiation becomes stronger 
under reduced gravity conditions. Therefore, detailed soot formation and radiation 
must be included in the flame structure analysis. A study of sooting turbulent 
diffusion flames under reduced-gravity conditions will not only provide necessary 
information for such practical issues as spacecraft fire safety, but also develop better 
understanding of fundamentals for diffusion combustion. 

Experimentally, full-field laser light transmission and thermophoretic soot particle 
sampling techniques will be used to measure flame soot particulate size and number 
density. Flame temperature will be measured using full-field two-color pyrometry 
and fine-gage thermocouples. These experiments are conducted in a drop tower. 

On modeling, the focus is on complete coupling of flame structure, soot formation, 
and radiation. The conserved scalar approach based on Favre-averaged governing 
equations, k-e-g turbulence model, and an assumed probability density function is 
used to predict flow field and gaseous species mole fractions profiles (refs. 3 and 4). A 
soot formation model developed by Syed et al. (refs. 5 and 6) will be modified to 
predict soot particle volume fraction and number density. The energy equation is 
incorporated to provide a full coupling between flame structure and radiation. The 
radiative heat flux is calculated from the radiative transfer equation (RTE) for a 
finite axisymmetric cylindrical enclosure, which will be solved by the P 3 spherical 
harmonics approximation (ref. 7). 

In this paper, a summary of the work to date and of future plans is reported. 
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Experimen tal Results 


Experiments are conducted using a 2.2-second drop rig. The rig is now instrumented 
for thermophoretic soot particle sampling and full-field laser light transmission 
measurements. Thermocouple and two-color pyrometry measurement are planned. 
Details for experiments are reported in another paper in this workshop (ref. 8). 


Results were collected for laminar flames first to validate our setup and procedure. 
The nozzle has a diameter of 1.7 mm. Propane and ethylene at 1.0 and 1.5cc/sec are 
tested. Micrographs of soot aggregates in 1.5cc/sec propane-air diffusion flames are 
shown in Fig. 1. The significant difference in aggregate size (or number of primary 
particles per aggregate) between normal and reduced gravity conditions is noted, and 
this has been observed for all fuel and flow rate combinations. 

Fig. 2 shows the mean diameter of primary soot particles as a function of height 
above nozzle exit. The size distribution, deduced from more than 100 particles per 
sample, agrees well with the normal distribution and 95% of the population falls in 
±10 nm range. Incipient soot particles are larger under 0-g. This seems to support the 
argument that incipient soot starts as liquid droplets, whose dynamics is affected by 
the relative dominance of surface tension. Under 0-g, the oxidation stage seems to be 
much weaker for the 1.5 cc/sec case, so flames may smoke (i.e., release soot particles). 
It is therefore concluded that rates of nucleation, growth, coagulation, and oxidation 
are very different between 1-g and 0-g. 0-g data are needed for determining these rate 
constants in models. Larger aggregate and primary particle sizes mean that radiation 
heat transfer is even more significant under 0-g (ref. 9). 

The full-field laser light transmission experiment has been tested for a 1-g flame. 
Soot volume fraction for a section as tall as 4.5 cm was measured instantaneously 
using a CCD camera. Results agree well with point-by-point measurements. 

Modeling of Turbulent let Diffusion Flame 


Following the conserved scalar approach, the structure of an axisymmetric turbulent 
diffusion flame is modeled using Fa vre- averaged boundary layer flow equations for 
conservation of mass, momentum, and mixture fraction described by Bilger (ref. 10), 
and a k-e-g turbulence model proposed by Lockwood and Naguib (ref. 11). The 
governing equations can be written in a general form as (refs. 12) 


d 1 d 1 d 

—(pu<p) + -—(rpv<t>) = -— 
dx r dr r dr 


d<f> 

^Tr 


+ S. 


( 1 ) 


where <j> = 1 (continuity), u (axial velocity), / (mixture fraction), k (turbulence kinetic 
energy), £ (turbulence dissipation), or g (mixture fraction variance). Favre-averaged 
(mass weighted) quantities are defined as <j>=p<(>/p, whereas an overbar represents 
conventional time-averaging. Details for (effective viscosity) and 5, (source 
term) can be found in the references. Buoyancy effects will only be considered in the 
mean flow equation, neglecting buoyancy-turbulence interactions. 
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Assuming all instantaneous scalar properties are functions of the mixture fraction 
only, termed the state relationships, the Favre-averaged mean and variance of scalar 
properties can be determined from an assumed probability density function (pdf). 
The mixture fraction is described as the mass fraction of fuel atoms anywhere in the 
flame and by virtue of its definition it is invariant with chemical reaction (so it is a 
conserved scalar). This allows the fluid mechanics to be decoupled from chemical 
reactions, and it is a good approximation if the rate of chemical reaction is limited by 
turbulent mixing and not by chemical kinetics (i.e. fast chemistry at equilibrium) 
and the diffusivities of all species and heat are equal. 

The system represented by Eq. (1) is solved using a block-tridiagonal code written by 
Chen (ref. 13). Major gaseous species concentrations are then calculated based on the 
state relationships constructed from thermodynamic equilibrium (adiabatic flame) 
calculations using STANJAN (ref. 14) and a /J-pdf. This has been found accurate for 
nonluminous and fuel-lean regions of luminous flames, but fails for rich regions of 
luminous flames. For better accuracy, the laminar flamelet approach (refs. 15 and 
16), with the state relationships constructed from laminar flame data, may be used. 

Modeling of Soot Formation and Oxidation 

Compared to gaseous specious, soot inception is a much slower reaction, and soot 
particles have very different diffusivities. In addition, soot particles are subjected to 
turbulent mixing, thermophoretic forces, and strong radiation heat transfer effects. 
The structure analysis outlined above has been extended to model soot formation 
and oxidation. Magnussen and Hjertager (ref. 17) introduced rate equation models 
for both nucleation (based on number density) and surface-growth (based on mass 
concentration) stages. They used an energy equation in the form of Eq. (1), with the 
source term for radiation, which is calculated using a two-flux model. Gore and 
Faeth (ref. 4) constructed state relationships for soot volume fractions from laminar 
flame data, and for temperature from equilibrium combustion calculation with a 
fixed fraction of chemical energy release lost by radiation. Kent and Honnery (ref. 18) 
found that the state-relationship approach for soot volume fractions may work in 
the top oxidation portion of the flame, where turbulent mixing dominates, but not 
the lower inception portion, where chemical kinetics dominates. They calculated 

the temperature from an energy equation with a radiation term of [~ecr(T* -T*)Jl]. 

It is believed that the rate-equation model approach is more accurate for predicting 
soot formation and oxidation. Kennedy et al. (refs. 19 and 20) applied Eq. (1) for soot 
volume fraction in laminar diffusion flames. The soot volume fraction source term 
was replaced by rate equation models derived for nucleation, growth, and oxidation 
with an assumed average number density. An energy equation similar to that in 
Kent and Honnery (ref. 18) was used. The model predicts peak soot volume fractions 
fairly accurately, but not their radial distributions. We adopt a two-equation model, 
derived by Moss et al. (ref. 6) and Syed et al. (ref. 5), for soot number density ( N ) and 
volume fraction (f v ). The model agrees well with measured flat (Wolfhard-Parker) 
diffusion flame data. The model has been modified (ref. 21) to couple with reduced 
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chemical reaction mechanisms, and successfully predicted radiation heat transfer 
from a turbulent reacting jet in a cross-wind (ref. 22). 


Since the model applies Eq. (1) to describe the transport of soot particles in terms of 
Favre-averaged number density and volume fraction, it can be easily incorporated 
into the reacting turbulent free-shear flow code. Rate equation models, representing 
nucleation, surfaced growth, coagulation, and oxidation, are used in source terms as 
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where n o = 6xl0 26 is Avogadro's number and p, is the mass density of solid carbon 

(typically 1.8- 2.0 g/cm 3 ). The rate constants a, /?, y, and % ar e modeled in terms of 
mixture density p , activation temperatures T a and T r , and fuel mole fraction X c , as 


a = C a p 2 T U2 X c exp(-T a /T) , (4) 

P = Cp T m , (5) 

y = C r pT m X c cxp(-T r /T), (6) 

X = C x T m X ox txp(-T z /T), (7) 

wherein coefficients and temperatures are determined from experimental data. In 
Eq. (7), X„ — Xpff if / > 0.064 (fuel-rich), and X 0I = 0.045X Oj if / < 0.064 (fuel-lean). The 
effective viscosities are given in the form of 


_ p, _ p Cu k 

a, ■ 

where G t are determined from experimental data and is the same as that for u . 
Flame Radiation Heat Transfer .Calculation 


( 8 ) 


Radiation is fully coupled with flame structure analysis through an energy equation 
given in the form of Eq. (1) with <f> = H (total enthalpy) and S $ ~-^ Q rai - Turbulence- 
radiation interactions are neglected. Gore et al. (ref. 23), using a multi-ray method to 
calculate the radiative heat flux q rad , showed that the coupled analysis predicts more 
accurately than the uncoupled one. However, their work neglected the scattering by 
soot particles, but accounted for the effects of turbulence-radiation interactions. 

The radiative transfer equation (RTE) for axisymmetric finite cylindrical enclosures 
such as for jet diffusion flames can be expressed as (ref. 7) 


if pjL-R d 4 ,,i 

fiy* dr r d<f> ^ dz 
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Here, the subscript A denoting spectral quantities has been left off, / is intensity (with 
subscript b for blackbodv), [5 is the extinction coefficient, co is the single scattering 
albedo, £,77, and fj. are direction cosines, and <P is the phase function. Both (5 and co are 
functions of position ( r , z) for axisymmetric diffusion flames. 

Menguc and Viskanta (ref. 7) derived a solution using the P 3 spherical harmonics 
approximation with the delta-Eddington approximation for the phase function. For 
soot agglomerates, the phase function can be better approximated by a third order 
Legendre polynomial series (ref. 9) given as 

Q(d,<p,d',<p') = l + ^a.P.fcosyO , (10) 

where y/is the scattering angle and cos y/ = 7777' + pfi', and coefficients a H are to be 

determined. Applying the spherical harmonics approximation, the model equations 
are obtained by employing the following integrations 

[Equation (9)]K."’sin0d9 dp , (11) 

for n- 0,1, N and m =-n, -n+ 1, — n, with N being the order of the P N approximation. 
The spherical harmonics Y" (superscript * denotes the complex conjugate) can be 
replaced with multiples of direction cosines. Define the moments of intensity as 

/ 0 (r,z) = J Q J q I(r,Q,<p,z)sinQd9 dp , (12a) 

t t )Kr,d,<p,z)sindd9dp , (12b) 

where each of the direction cosines l t , and 4 is £,*7/ or P- ■ Physically, the zeroth 
moment is the total incident radiation, whereas the first moments are radiative heat 
fluxes along the coordinate axis. The divergence of radiative heat flux vector, which 
appears in the energy equation, can then be calculated from 

V ■ 4n* = Pd - <0){4nI h [T(r,z)]- I 0 (r,z)} = S(r,z) . (13) 

After some manipulations, results obtained from Eq. (11) can be combined into four 
coupled elliptical partial differential equations for 7 0 , /n, / 33 , and / 13 . The model 
equation for is given as 

+ r„) + i(4B, + 7B, )r, + + 14S,) + 1^B 2 T 1 

where B's are functions of a n and co, r = (5r 0/ Ai is a function of B, z, I n , / 33/ and / 13 , and 
T is the derivative operator. These model equations, together with 16 Marshak's 
boundary conditions derived for diffusely emitting and reflecting opaque boundary, 
are solved numerically using a solver called ELLPACK (ref. 24). 


, , = -fr I | + A. H4) 
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To incorporate radiation calculations into the overall flame structure analysis, there 
are two major difficulties to be overcome. First is that equations for the moment of 
intensity such as Eq. (14) are elliptical, whereas equations represented by Eq. (1) are 
parabolic. It seems impossible to solve both systems simultaneously- An alternative, 
which we are testing now, is to iterate between these two systems until the results 
numerically converged. The second difficulty is on modeling the spectral radiative 
properties for integrating the RTE over the spectral range. We will first follow the 
wide-band approach used by Song and Viskanta (ref. 25). 


Preliminary Result s and Discussion 


We first tested the flame structure and soot formation /oxidation predictions against 
results in Kent and Honnery (ref. 18) for an ethylene-air turbulent diffusion flame 
(nozzle diameter D =3 mm, exit velocity U 0 = 52 m/sec, Re 0 =9615). Good agreement 
was found for mixture fractions. Fig. 3 shows the comparison for the temperature 
along the centerline and radially atx/x m = 0.4,0.7Cr m = 115D). The dashed curve is from 
adiabatic combustion calculation, and the solid curve from a correction (ref. 21) 


T = T 



(15) 


where T ad is the adiabatic temperature, T mai is the maximum adiabatic temperature, 
and 0.09 <b< 0.15. Even Eq. (15) still over-predicts at the flame tip and outer edge. 
This demonstrates the importance of including the energy equation and detailed 
radiation heat transfer calculations. 


Fig. 4 shows that the soot formation/oxidation model predicts soot volume fraction 
accurately. The over-prediction near the outer edge is caused by that of temperature. 
Fig. 5 shows the temperature and soot volume fraction contours, calculated using 
the same rate coefficients, for one of our test conditions {Reo =536, Re L = / 400) . The 
temperature contour shows a much taller and wider 0-g flame, as expected. The 0-g 
flame soot volume fraction contour is obviously unreasonable, since different rate 
constants are needed as concluded from laminar flame data. A reasonable prediction 
is possible by decreasing the level of oxidation for 0-g, which agrees with laminar 
flame results. 0-g measurements are needed for determining the rate constants. 

The implementation of a iterative proce ss betw een the RTE and the flam e s tructure 
solvers has not been completed. Some numerical problems have been experienced, 
which are believed to be caused by zero /3 in part of the cylindrical enclosure where 
there is no flame. Fig. 6 shows the results of a forward calculation of the terms in the 
energy equation. Referring to Eq. (1) with <f> = H and S # = — V ■ those two terms on 
the left are convection terms, and the first term on the right is conduction with the 
last term being radiation. The first three terms are calculated using total enthalpy 
from the adiabatic combustion calculation for 1-g. The last term is calculated using 
the temperature from Fig. 5 with an assumed uniform soot volume fraction of 10 8 . 
It is clear that all these terms are of about the same magnitude, and the radiation 
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term depends strongly on location. This confirms the importance of incorporating 
detailed radiation calculation into the flame structure analysis. 

The modeling can be improved by accounting for interactions between turbulence 
and radiation (refs. 12 and 25), and between turbulence and buoyancy (ref. 22). The 
latter should be an important factor between 1-g and 0-g conditions. 
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(b) 0-g flame, 80 mm HAB 1 JJ.DTI 


Figure 1. Micrographs of soot aggregates in 1.5 cc/sec propane-air laminar diffusion 
flames under (a) normal, and (b) reduced gravity conditions. (1.7 mm nozzle dia.) 
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Figure 2. Mean diameter of primary soot particles in laminar diffusion flames. 
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Figure 3. Comparisons of predicted and measured turbulent flame temperatures 
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Figure 6. Estimated values of various terms in the energy equation with radiation. 

132 



COMMENTS 


Question (J.P. Gore, Purdue University): 

(1) What is the magnitude of turbulence-radiation interactions, and what is the justification for neglecting 
this effect in preference to scattering? Scattering is more important in optical diagnostics in the visible 
rather than in the thermal radiation term in energy equations. 

(2) The effects of joint probability distribution functions between mixture fraction and enthalpy arc 
important. Do you plan to model these? 

(3) Are the soot volume fraction data from Kent and Honncry (1987) path-integrated equivalent 
measurements? If so, how is the comparison with predictions made? 

(4) Chemical equilibrium calculations arc not in agreement with measurements on the fuel-rich side. Should 
the model be improved in this aspect? 

(5) In your verbal response you mentioned that there arc other uncertainties such as refractive indices of 
soot. However, work reported in the literature has shown that items (1), (2) and (4) above have a very 
strong (even order of magnitude) influence on radiation heat loss. In my opinion, these need to be 
addressed with priority. 

Answer: (1) There is no attempt at this point to model or estimate the magnitude of turbulence-radiation 
interaction, or to justify that it is negligible. The turbulence-radiation interaction is not included due to its 
complexity and the lack of definitive information about its effects. On the other hand, the scattering effects 
have been shown to be significant (Ref. 9) and can be dealt with in our radiation heal transfer calculations. 

(2) We agree that the effects of joint pdf between mixture fraction and enthalpy are important and should be 
included in the model. However, this will not be our focus. Instead, we plan to follow works in existence, 
such as those by the commentator (Ref. 23), and those developed in the future. 

(3) We believe that soot volume fraction data in Kent & Honnerv (Ref. 18) are local values, instead of path- 
integrated. The last sentence in the second paragraph in Experimental Procedure stales that horizontal 
traversing and Abel inversion were applied to extract local extinction data. 

(4) Yes. We will improve this aspect of the model by using t he laminar (lamelel model, as indicated in the 
paper. 

(5) One may argue that soot refractive index data are fairly accurate in the visible wavelengths, even though 
we are not perfectly convinced of that. However, it is fair to say that data in the infrared arc highly 
questionable. Since there is no directly related analysis on the effects of this uncertainly, we plan to 
investigate this aspect more quantitatively. 


Question: (Ivan Cation, UCLA): In transport processes, one is usually interested in number density and a 
mean diameter. Your electron micrographs show shapes as far from spherical, and one has to wonder what 
the "mean diameter" represents. 

Answer : Wc consider soot particulates as aggregates of near-spherical primary particles. The shapes, as 
referred to in the question, are shapes of aggregates. The mean diameters, as referred to in the paper, are 
mean diameters of the near-spherical primary particles. 


133 








